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Abstract 

A calculational scheme is developed to evaluate chiral corrections to prop- 
erties of composite baryons with composite pions. The composite baryons 
and pions are bound states derived from a microscopic chiral quark model. 
The model is amenable to standard many-body techniques such as the BCS 
and RPA formalisms. An effective chiral model involving only hadronic de- 
grees of freedom is derived from the macroscopic quark model by projection 
onto hadron states. Chiral loops are calculated using the effective hadronic 
Hamiltonian. A simple microscopic confining interaction is used to illustrate 
the derivation of the pion-nucleon form factor and the calculation of pionic 
self-energy corrections to the nucleon and A(1232) masses. 
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I. INTRODUCTION 



The incorporation of chiral symmetry in quark models is an important issue in hadronic 
physics. The subject dates back to the early works [0-0 aimed at restoring chiral symmetry 
to the MIT bag model . The early attempts were based on coupling elementary pion fields 
directly to quarks. A great variety of chiral quark-pion models have been constructed since 
then and the subject continues to be of interest in the recent literature Despite the 

long history, there are many important open questions in this field. In the present paper, 
we are concerned with one of such questions, namely the coupling of the pion as a quark- 
antiquark bound state to the baryons. Starting from a model chiral quark Hamiltonian, we 
construct an effective low-energy chiral pion-baryon Hamiltonian appropriate for calculating 
chiral loop corrections to hadron properties. The composite pion and baryon states are 
determined by the same underlying quark chiral dynamics. 

The model we use belongs to a class of quark models inspired in the Coulomb gauge QCD 
Hamiltonian [|ll|] and generalizes the Nambu-Jona-Lasinio model Jl2| to include confinement 
and asymptotic freedom. This class of models is amenable to standard many-body techniques 
such as the BCS formalism of superconductivity. The initial studies within these models 
were aimed at studying the interplay between confinement and dynamical chiral symmetry 
breaking (D%SB), and concentrated on critical couplings |TB[ for D%SB and light-meson 
spectroscopy fl4}| . The model has been extended to study the pion beyond BCS level and 



meson resonant decays in the context of a generalized resonating group method ||15|| . Since 
the model is formulated on the basis of a Hamiltonian, it provides a natural way to study 
finite temperature and chemical potential quark matter |16|]. The model and the many-body 
techniques to solve it make direct contact with first-principle developments such as nonper- 
turbative renormalization-group treatments of the QCD Hamiltonian Jl7| and Hamiltonian 
lattice QCD £§]. 

One important development of the model, in the context of the present paper, was 



its extension in Ref. [19| to baryon structure. In Ref. [19] a variational calculation was 
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implemented for the masses of baryons and it was shown that a sizeable A(1232)-iV mass 
difference is obtained from the same underlying hyperfine interaction that gives a reasonable 
value for the ir-p mass difference. This hyperfine interaction, along with other spin-dependent 
interactions like tensor and spin-orbit, stem from Bogoliubov-Valatin rotated spinors that 
depend on the "chiral angle" . The chiral angle gives the extent of the chiral condensation in 
the vacuum and determines the chiral condensate. The very same variational wave function 
was used later for studying S-wave kaon-nucleon |20| scattering and the repulsive core of 



the nucleon-nucleon force ||21|| . Both calculations obtain S'-wave phase-shifts that compare 



reasonably well with experimental data. A remarkable feature of all these results for the 
low-lying spectrum of mesons and baryons and S'-wave scattering phase-shifts is that they 
are obtained with a single free parameter, the strength of the confining potential. 

In the present paper we go one step forward in the development of the model and set 
up a calculational scheme to treat chiral corrections in hadron spectroscopy. In a recent 
publication |[22|| , two of us have calculated the pion-nucleon coupling constant in this model 



and obtained reasonable agreement with its experimental value. Here, we are interested in 
developing a scheme for calculating chiral corrections to hadron properties. We study the 
requirements to obtain in the context of the model the correct leading nonanalytic behavior 
(LNA) of chiral loops. Our scheme follows the standard practice f| f23| of constructing an 



effective baryon-pion Hamiltonian by projecting the quark Hamiltonian onto a Fock-space 
basis of single composite hadronic states. Chiral loop corrections are then calculated with 
the effective Hamiltonian in time-ordered perturbation theory. The difference here is that 
while in the previous works the pion is an elementary particle, in our approach the pions 
are composites described by a Salpeter amplitude. 

A difficulty appears in the implementation of the projection of the microscopic quark 
Hamiltonian onto the composite hadron states, which is not present when the pion is treated 
as an elementary particle. The difficulty is related to the two-component nature of the 
Salpeter amplitude of the pion. The two components correspond to positive and negative 
energies (forward and backward moving, in the language of of time-ordered perturbation 



theory), they are 2x2 matrices in spin space and are called energy-spin (E-spin) wave 
functions. For the pion, the negative-energy component is as important as the positive- 
energy one - in the chiral limit they are equal - because of the Goldstone-boson nature of 
the pion. Because of this, the Fock-space representation of the pion state is not simple. 
We overcome the difficulty by rephrasing the formalism of the Salpeter equation in terms 
of the RPA (random-phase approximation) equations of many-body theory. The single- 
pion state is obtained in terms of a creation operator acting on the RPA vacuum. The 
pion creation operator is a linear combination of creation and annihilation operators of 
pairs of quark- ant iquark operators; the positive-energy Salpeter component comes with the 
creation operator of the quark-antiquark pair, and the negative-energy one comes with 
the annihilation operator of the quark-antiquark pair. In this way, the projection of the 
microscopic quark Hamiltonian onto the single hadron states becomes feasible and simple. 

The paper is organized as follows. In the next Section we review the basic equations of 
the model. We show the relationship of the formalisms of the Salpeter equation and of the 
many-body technique of the RPA. In Section [ITJ we derive the pion-baryon vertex function 
in terms of the bound-state Salpeter amplitudes for the pion and the baryon. We obtain an 
expression that is valid for a general microscopic quark interaction, not restricted to a specific 
form of the potential. Given the pion-baryon vertex function, we derive the expression for 
the baryon self-energy correction in Section |iy]. Although the derivation of the expression 
for the self-energy is well-known in the literature, we repeat it here to make the paper easier 
to read. In Section [V| we obtain numerical results for the pion-baryon form factor and 
coupling constants. Numerical results and the discussion of the LNA contributions to the 
baryon masses are presented in Section [VI]. Section [VIl] presents our conclusions and future 



directions. 



II. THE MODEL 



The Hamiltonian of the model is of the general form 
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H = H + H r , (1) 

where H Q is the Dirac Hamiltonian 

H = J dxip^(x) (-ict-V + (3m q )ip(x), (2) 

with ip(x) the Dirac field operator, and Hj a chirally symmetric four-fermion interaction 

H I= l -jdxj dy^(x)T a T^(x)V T (x-y)^(y)T a r^(y). (3) 

Here, T a = 1/2 A a , a = 1, . . . , 8 are the generators of the color SU(3) group, T is one, or 
a combination of Dirac matrices, and Vr contains a confining interaction and other spin- 
dependent interactions. One example of Vr will be presented in section |V|, when we make a 
numerical application of the formalism. 

Once the model Hamiltonian is specified, the next step consists in constructing an explicit 
but approximate vacuum state of the Hamiltonian in the form of a pairing ansatz. This is 
most easily implemented in the form of a Bogoliubov-Valatin transformation (BVT). The 
transformation depends on a pairing function, or chiral angle (p that determines the strength 
of the pairing in the vacuum. The quark field operator is expanded as 

4>{*) = J J0ff 2 [«.(«) Kq) + v s (q) <ft{-q)] , (4) 

where the quark and antiquark annihilation operators b and d annihilate the paired vacuum, 
or BCS state, |0bcs)- Here the spinors u s (q) and v s (q) depend upon the chiral angle (p as 

u s (v) = -j= {t 1 +sin^(g)]i + [1 - sin^(g)]^a ■ q] u° s , (5) 

v s (q) = -^{[l+sin^(g)]^ - [1 - sin (p(q))^ a ■ q] v° s , (6) 

where u° s and v® are the spinor eigenvectors of Dirac matrix 7 = (3 with eigenvalues ±1, 
respectively. 

The chiral angle can be determined from the minimization of the vacuum energy density, 

Evac f dq r 1 
= / - — -77 lr ct ■ q A (q) 

V J (2tt) 3 L y Vy; J 

+ lfwy^ {q - 0* [rrA Hm - W ] , (7) 



where Tr is the trace over color, flavor, and Dirac indices, and 



A + (<7) = E«-(«)«!(«) > A ^<?) = T,Vs(qW s ((q). 

s s 

The minimization of the vacuum energy leads to the gap equation, 

A(q) cos (p(q) — B(q) sin (p(q) = 0, 

with 



(8) 



(9) 



(10) 
(11) 



A(q) = m + yWv T (q- q') Tr [(3T a T (A+(,/) - K~(q)) T a T 
B{q) = q + \ J ^ V T {q ~ q')Tr [a • qrTT (A+(q') - A"( 9 )) TT 

where Vr(<?) is the Fourier transform of Vr(aj), 

\7 r (q) = y daje^^te)- (12) 

The pion bound-state equation is given by the field-theoretic Salpeter equation. The wave 
function has two components, + and <p~, the positive- and negative-energy components. 
Each of the 0's is a 2 x 2 matrix in spin space, 0+ S2 and 0~ S2 . For this reason, the 0's 
are also called energy-spin (E-spin) wave functions. The 0's satisfy the following coupled 
integral equations: 



[M(q) - E{q + ) - E(q^+(q) = u\q + ) K^k, q) v(q_), 
[M(q) + E{q + ) + £(g_)]0 fc T (q) = v\q + ) K^k, q) u(q_), 

where q± = q ± k/2, and the kernel K^k, q) is given by 

K,(k, q) = J -^Mq - q') T a r [u( q ' + ) <t>t{ q >) v\k'_ 

+ v{k' + )<P k T {q>)u\q'_)]T«T. 



(13) 
(14) 



(15) 



Here, the super-script T on means spin-transpose of <f> , (f> s * S2 = 4 ) S2 ,s 1 - The amplitudes 
are normalized as 
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/ Wi^Uq) - KWtwla)] = - *0 ■ (16) 

In the language of many-body theory, these equations can be identified with the RPA 



equations [24]. In the RPA formulation, one writes for the one-pion state (suppressing 



isospin quantum numbers) 

|7r(fc)) = Mt(fc)|0 RPA ), (17) 

where |0rpa) is the RPA vacuum, which contains correlations beyond the BCS pairing 
vacuum |0bcs), an d M\{k) is the pion creation operator 

M\{k) = J j0^j- 2 6t(g + ) <#{qj) - <fc(q) b(q + )d(q_)\ , (18) 



where the q± were defined above. In this formulation, Eqs. Q13D and (|14|) above are obtained 
from the RPA equation of motion: 

(7r(fc)|[#,il4]|0 R p A ) = (E w (k) - E vac ) (ir(k)\M^\0 RPA ). (19) 

The normalization is such that 

(7r(fc)k(fc')> = (OR,PA|[M^(fc),Mt(fc , )]|0 R p A ) = 6(k-k'). (20) 

The verification of D^SB consists in finding nontrivial solutions to the gap equation, 
Eq. (|9|), and the existence of solutions for the pion wave function. Refs. [11) and |I5| show 



that for a confining force, there is always a nontrivial solution for the gap and pion Salpeter 
equations. Moreover, good numerical values are obtained for the chiral parameters such 
as the pion decay constant and the chiral condensate when an appropriate spin-dependent 
potential is used |25[| . 

The inclusion of RPA correlations beyond BCS pairing was shown in Ref |24| to have 
a dramatic effect on the mass spectrum of the pseudo scalar mesons (n and if), while it 
has almost no effect on the mass spectrum of vector mesons (p and u). One can trace this 
effect to the fact that the pseudo scalar mesons have a sizeable "negative energy" component 



wave function, while the vector mesons have a very small negative energy component |T5 . 
For baryons (such as nucleon and A), since they do not have a sizeable negative energy 
component one expects that they can be reliably obtained from the BCS vacuum. We 



write therefore for the one-baryon state 

|£(°)(p)) = £(°)t(p)|0 BCS >, (21) 

where the baryon creation operator B^(p) is given by 

B<°)t(p) = J dq x dq x dq x S(p - q x - q 2 - q 3 ) ^ P (<Ziq 2 <? 3 ) e ClC2C3 x ^ s ^ 

x&L 1 / 1 (9i) & L/ 2 (92)&L/a(93). (22) 

Here e ClC2C3 is the Levi-Civita tensor, which guarantees that the baryon is a color singlet and 
-^/i/2/3sis2S3 are tjjg spin-isospin coefficients. The wave function ^! P (qiq2qz) is determined 



variationally [19|. The index (0) on the baryon operators indicates a bare baryon, i.e. a 
baryon without pion cloud corrections. 

The important fact to notice here is that the baryon wave function depends on the chiral 
angle <p and as such, spin splittings and other properties are determined by the same physics 
that determines vacuum properties. The pion-baryon vertex, that we will discuss in the 
next section, will therefore depend on the chiral angle not only because of the pion, but also 
because of the baryon wave function. Numerical results for the masses of the the nucleon 



and A(1232) have been obtained previously ]19| and are of the right order of magnitude as 



compared with experimental values. Of course, fine-tuning with different spin-dependent 
interactions can improve the numerical values of the calculated quantities. 

III. PION-BARYON VERTEX FUNCTION 

We obtain an effective baryon-pion Hamiltonian by projecting the model quark Hamil- 
tonian onto the one-pion and one-baryon states, Eqs. ([T7|) and (EH|). We use a shorthand 
notation. For the bare baryons, i.e. baryons without pionic corrections, we use the indices 



S 



a, /?,••• to indicate all the quantum numbers necessary to specify the baryon state, such as 
spin, flavor and center-of-mass momentum. For the pion, we use j, k, ■ ■ ■ to specify all the 
quantum numbers of the pion state. With this notation, the effective Hamiltonian can be 
obtained as 

H = £ \a)(a\H\(3)((3\ + £ \j)(j\H\k)(k\ 

af3 jk 

+ E(i«)b'}(^il^il^)^l + \(3}(f3\Hi\j^)(a\(j\) . (23) 

jaf3 

This leads to an effective Hamiltonian that can be written as the sum of the single-baryon 
and single-pion contributions, and the pion-baryon vertex: 

H = H + W, (24) 

where Ho = Hb + H n contains the single-baryon and single-pion contributions 

H B = J2 B<£\ H« = £ E d M] M 3 , (25) 

a j 

and W is the pion-baryon vertex 

W = J2Wi,B^B^M ]+ h.c. (26) 

jot/3 

Here, and (Mj and Mj) are the baryon (pion) creation and annihilation operators, 
discussed in the previous section. Note that we have assumed that states with different 
quantum numbers are orthogonal. Note also that in writing the state \ j, a) we have implicitly 
assumed that the negative-energy component of the baryon is negligible and the baryon 
creation operator acting on the RPA vacuum has the same effect as acting on the BCS 
vacuum. 

We note that the projection of the microscopic quark Hamiltonian to an effective hadronic 
Hamiltonian can be obtained in a systematic and controlled way using a mapping proce- 



dure ||26|| . We do not follow such a procedure here because we are mainly concerned with 
tree-level pion-baryon coupling and the projection we are using is enough to obtain the de- 
sired effective coupling. For processes that involve quark-exchange, such as baryon-meson or 

9 



baryon-baryon interactions, a mapping procedure would be useful. In a future publication, 
we intend to address such processes in the context of the present model. 
The single-hadron Hamiltonians Hq and give 

Ho\B^) = Ef\B^) , ff» = , (27) 

with E$ = Mg and = in the rest-frame. The vertex W gives the coupling of 
the pion to the baryon and, as we will show later, leads to loop corrections to the baryon 
self-energy. 

The pion-nucleon vertex function can be written generally as 

W = j2(w+ + Wr) , (28) 

4 = 1 

where the W^s are of the general form (for simplicity we suppress the color and spin-flavor 
wave functions in the following) 

Wt(p,p';k) = J V r (q - q') ^(qW^'s) F, 0k] ^Aq^), (29) 

where the VV^'s involve the Dirac spinors and the pion wave functions. In Figure 1 we present 
a pictorial representation of the different contributions to the vertex function. Explicitly, 
the 's are given by 

W+ = [u t (gi)T a ru(g 4 )] <PUp±) \v\-p A )T a Yu{ qi )} , (30) 

W 2 + = <t>t(pA) [v^-p,)T a Tu( qi )} \u\q' 2 )T a Tu{q 2 )} , (31) 

W 3 + = W+ , (32) 

Wf = [u^q'^Tvi-q,)] fc r (p 4 ) [u^p A )T a Yu{ qi )\ , (33) 

W 2 = \u\q' 2 )T a Tu{q 2 )} [u^q'^Tv^q,)] fc r (p 4 ) , (34) 

VV 3 _ = W 2 . (35) 

In these formulas, the quark momenta in the initial (final) nucleon q l5 q 2 , q% (q[, q' 2 , q' 3 ) and 
the momenta of the quark and antiquark in the pion, and q^, are expressed in terms of 
the loop momenta q, q', q" by momentum conservation (see Fig. [l]) . 
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Once the effective baryon-pion Hamiltonian is obtained, one can calculate the pionic 
corrections to baryon properties as in the CBM and in the traditional Chew-Low model. 
This will be done in the next section. 

Before leaving this section, we recall that the use of the Breit frame is essential in 
calculations of form factors (vertex functions) in static models [57], like the present one. 



This is true for composite models for which approximate solutions that maintain relativistic 
covariance are very difficult to implement. This was the case for all old, static source, 
pion-nucleon models like the Chew-Low model ||28|| . In particular, as explained in Ref. |27| . 



electromagnetic gauge invariance is respected in this frame. Therefore, in calculating loop 
corrections to baryon properties, we employ the Breit-frame vertex functions. In the Breit 
frame, we denote the incoming pion and nucleon momenta by p and —p/2, respectively, and 
the outgoing nucleon momentum by p/2. In this frame, the internal momenta of quarks and 
antiquarks q' 1} q' 2 , ■ ■ ■ are given in terms of the loop momenta q, q' and q" as 

Ql = p/2 + q + q" q[ = p/2 + q' + q" 

q 2 = q' 2 = -q' + q" q 3 = q' 3 = -2q" 

P*=p/2-q-q" q±=p/2-q-q\ (36) 
for the vertex W + and 

Ql = -p/2 + q' + q" q\ = p/2 + q' + q" 

q 2 = q ' 2 = -q> + q " q 3 = q' 3 = ~2q" 

p 4 = -p/2 + q + q" q i = - p /2-q-q\ (37) 

for the vertex W~. In following equations, we also denote the vertex function as 
W(-p/2,p/2;p) = W{p). 

IV. SELF-ENERGY CORRECTION TO BARYON MASSES 



For completeness we review the derivation of the expression for the self-energy correction 
from the effective baryon-pion Hamiltonian of Eq. (|4]) in the "one-pion-in-the air" approxi- 
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mation p8|,|6| . The baryon self-energy is defined as the difference of bare- and dressed-baryon 
energies: 



T,(E B ) = Eb — Eg\ 



The physical baryon mass Mb is given by the (in general, nonlinear) equation 



M B = il4 0) + S(M B ) 



(38) 



(39) 



where Mg is the bare baryon mass (i.e. without pionic corrections) and H(Eb) is the 
self-energy function. 

Let | B) denote the physical baryon state, and \B ) the "bare" undressed state. Let Z^ 
be the probability of finding \Bq) in \B). Then one can write 



B) = JZi\B )+A\B) 



where A is a projection operator that projects out the component \B ) from \B), 



A = 1-\B )(B \. 



(40) 



(41) 



We have that 



(B\W\B ) = (B\(H - H )\B ) = (E B - eP)(B\B ) 
= Jzf E{E B ) . 



%2 (Eb — E B 



(42) 



We can now express \B) in terms of \Bq) and the pion-baryon interaction Hamiltonian W 

as 



\B) = JZ? 



1 



E B -H ~ AW A 



W 



\Bo). 



(43) 



On the other hand, since 



(B\W\B ) = JZ^(B \W 



we have that the self-energy is given by 



Eb — Hq — AW A 



W\B ) 



(44) 
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E(E B ) = -±= (B\W\B ) 



= ^ W E B -H -AWA W ^- (45) 

This expression can be further approximated so as to avoid solving complicated integral 
equations for the self-energy. We can manipulate the expression for E to obtain (for details, 
see Ref. @): 

E(E B ) = (B \W 1 W\B ) , (46) 

&B — -"0 — 2->0{-&b) 

with 

S (E B ) = WA 1 - AW. (47) 
The approximation consists in absorbing £o (-£?#) into Hq such that 

H + i: (E B ) = H , (48) 

with 

H = J2 E « Bi " )] 4° } + E M) Mj, (49) 

a 3 

where E a and .E^ are the physical energies. Therefore, the baryon self-energy can be written 

as 

E{E B ) = (B \W 1 W\B ) . (50) 
Eb — B 

Finally, insertion of a sum over intermediate baryon-pion states in Eq. (|50D leads to 

E(E B ) = Y / (B \W\n)—^-—(n\W\B ). (51) 

n — &n 

The structure vertex-propagator- vertex W(E—H )~ 1 W in Eq. ([)(]) is an effective baryon- 
pion interaction. The main difference here with the hybrid approaches |1|-§| is that we do 
not have a point-like pion coupling to point-like quarks and antiquarks. The pion-baryon 
vertex arises through the "Z-graphs" in which the antiquark of the pion is annihilated with 
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a quark of the "initial" baryon and the quark of the pion appears in the "final" baryon. 
Therefore, the vertex function incorporates not only the extension of the baryons, but also 
the extension of the pion. 

We truncate the sum over the intermediate states in Eq. (|51"D to the lowest mass states, 
namely the nucleon and the A(1232). In this case, we obtain for the on-shell N and A(1232) 
self-energies the coupled set of equations 

dk 



X N (M N ) 
£ A (M A ) 



(2tt) 3 
dk 



(2 



7T 



\W NN (k)\ 2 W NA (q)W AN (k) 
M N - (M N + E n (k)) + M N - (M A + E n (k)) 
W AN (k)W NA (k) |^ AA (fc)| 2 



(52) 
(53) 



M A - (M N + E^k)) M A - (M A + E w (k))_ 

This is the final result for the pion loop correction for the nucleon and A(1232). 

One important consequence of projecting the microscopic quark interaction onto hadronic 
states is that the leading nonanalytic (LNA) contributions in the pion mass as predicted 
by chiral perturbation theory are correctly obtained |J. In particular, as we discuss in the 
next section, Eq. ( [52] ) leads to an LNA contribution to the nucleon mass as predicted by 
QCD [f2g], namely 

m - na = -t^t/^- (54) 

In models where the pion is treated as a point like particle, this result follows trivially 
from Eq. fl52"|). In the context of the present model, where the pion is not treated covariantly, 
such a result does not follow in general for an arbitrary interaction. The difficulty is related 



to the fact the the pion dispersion relation, E w = y k 2 + m 2 , is not obtained in general in a 
noncovariant model. In the CBM for example, the pion is point like and the normalization is 
correct from the very beginning. However, the microscopic quark interaction can be chosen 
such that the pion dispersion relation is correctly obtained |n,2"5|. These issues will be 
discussed in Section |VT[ 
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V. THE PION-NUCLEON AND PION-A(1232) FORM FACTORS 



Our aim is to obtain an estimate for the numerical values of the pionic self-energies. 
It happens that nature has produced a sort of low energy filter (chiral symmetry) for the 
details of strong interactions. Indeed it is remarkable that although intermediate theoretical 
concepts like gluon propagators, quark effective masses and so on, might vary (in fact they are 
not gauge invariant and hence they are not physical observables) , chiral symmetry contrives 
for the final physical results, e.g. hadronic masses and scattering lengths, to be largely 
insensitive to the above mentioned theoretical uncertainties. The pion mass furnishes the 
ultimate example: In the case of massless quarks, the pion mass is bound to be zero, 
regardless of the form of the effective quark interaction provided it supports the mechanism 
of spontaneous breakdown of chiral symmetry. The other example is provided by the ir — n 



scattering lengths |Kj which are equally independent of the form of the quark kernel |3T 
Furthermore it has become more and more evident through the accumulation of theoretical 
calculations on low energy hadronic phenomena, ranging from calculations on Euclidean 
space to instantaneous approximations and from harmonic kernels to linear confinement, 
that low energy hadronic phenomenology only seems to depend mildly on the details of the 
quark kernels used. To this extent, we will use for the quark-quark interaction a kernel of 
the form, 

r = 7 °, (55) 

V(k) = ^(2vr) 3 K 3 A fc( 5(k), (56) 

where K is a free parameter. This potential has been widely used in the context of chiral 
symmetry breaking because it allows a great deal of simple analytic calculations (which is 
not the case for the linear potential). The harmonic potential basically differs from the 
linear potential in domains of the baryon-pion-baryon overlap kernel which contribute little 
to the total geometrical overlap so that, at least for results proportional to these overlaps, 
they should not differ too much. 
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The momentum-dependent part of the Salpeter amplitude for the baryon, \I/ p (gi<Z2<?3) 
in Eq. fl22]), is taken to be of a Gaussian form 

,1, („ „ „ 1 e-'"" +A ''/ 2 °° P1-P2 , Pi + P2 - 2p 3 .... 

where a 2 B is the variational parameter and Mb{p) is the normalization. Notice that since 
the integrations of the quark momenta in the functions Wf" in Eq. ( p9|) are made through 
a Monte Carlo integration, the Gaussian ansatz is not essential and does not simplify our 
calculations, but we still use it to make contact with previous calculations. 



As in our previous calculation [22] for the pion-nucleon coupling constant, the Salpeter 
amplitudes 4>k(q) up to first order in k are given by 

4>k(Q) - Nik)' 1 [+siny?(q) + Ei(fc)/i(g) + i gi (q) k-(q x a)} X ,S color , (58) 
<l>h{Q) - Nik)' 1 [- sin (p(q) + E 1 {k)f 1 {q) -ig 1 {q)k'{q x a)] x^S^ , (59) 

where if is the chiral angle and E\(k) is the first-order correction to the pion energy. The 
normalization M{k) is proportional to Ei(k) and is given as 

M 2 {k) = m(k) I ^3 sin^(g) h{q) = E x a\ (60) 

The energy E\(k) is given in terms of the second derivatives of the diagonal components of 



the Salpeter kernel with respect to k and its explicit form is given in Eq. (24) of Ref. [22 
Note that the truncation up to first order in k of the Salpeter amplitude constitutes a 
reasonable approximation due to the fact that cm. momenta-dependent distortions of the 
pion and nucleon wave functions are geometrically damped because of the geometric overlap 
kernel integrations for the functions in Eq. ( |29| ) - see Ref. ||32j| . Explicit numerical 
solutions were obtained in Ref. E2] for the functions fi(q) and gi(q)- 



For completeness, we initially repeat the results of Ref. |22[ for the coupling constants 
fnNN an d fnNA- in Ref- 123, they were obtained as 



fnNN 5 Of s (p) A 

(T 'V'P = wS n cttvP, (61) 



3^ 2a 

' 2V2 OfM /n tysiP) 

y/3 2a 2a 



fnNA c, 



Sp, (62) 
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where the isospin matrix is omitted and 



o fs {v) 



J[dq] [a + + a~ + b + + b~) 



(63) 



where [dq] means integration over q, q' and q" (see Eq. fl29|) ) and the set of functions 



a + ,a ,b + ,b is given by, 

a + = <p + \ <p'(qi)tr ■ W out + 



v'(gi) + 



a = ^ out { (p'(q[)<T ■ V*i„ 4 
b + = 



cosy? fa) 
cosy? fa) 



9i 



1 - sin^>fa) J cos^(pi) 



2gJ 



Pi 



(T -qiX (qi x V^ ou t)si | * 4 

cr-qr; x (^ X V^ n )j</T, 

cosy^fa; 



?1 



1 -siny?fa) [ cosyfa)) 
o = z < 2. : (T ■ q 1 + 



2gi 

X ^out *m ' 



q'i 



cr-q 1 x{q l x q x ) 



(T-q 1 x{q l x q 1 



(64) 
(65) 

(66) 
(67) 



Here, ^i n> out stand for the baryon m and out Salpeter amplitudes and represent the 
pion Salpeter amplitudes. 

The baryon-pion coupling constants are obtained as the zero limit of the nucleon (or 
A) momentum, p — > 0, of the above overlap functions. For simplicity, we are defining the 
couplings at zero momentum, and not at the physical pion mass. In order to facilitate 



the integration, in Ref. |22| a Gaussian parameterization for the (cos (p(k))/k and (1 — 
sin (p(k))/k 2 was used. Here, since we need the vertex function for p ^ 0, we use a Monte 
Carlo integration to perform the multi dimensional integral that gives the overlap function 
and use the full numerical solution for the gap function (not the Gaussian parameterization). 
We first checked the correctness of our Monte Carlo integration with the result of Ref. |22[ 
for the special case of p = using the same Gaussian parameterization as was used there. 
This was done by calculating Of s (p) for p = (0,0, p z ) and finding the limit of O = Of s /p z 
when p z — > to obtain f n NN- 

As in Ref. we have used Kq = 247 MeV for the strength of the potential. The 
variational determination of a of the baryon amplitude, Eq. (]57|) , leads to = 1.2Kq. For 
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the A(1232), the result is not much different and therefore we use a at = OA- 
Introducing the quantities 

Fl = i7ir F ' = i7if- (68) 

we can summarize the couplings of the pion to the nucleon and A(1232) as follows: 

777 777 777 

f wNN = Fx O(0) — , fnNA = F 2 O(0) — , fnAA = F 3 O(0) — . (69) 
a a a 

For the value of K Q given above, we have m^/a = 3.47. The numerical values for the 
couplings are then 

U NN = 1.19, f nNA = 2.02, / wA a = 0.24. (70) 

The effect of the Gaussian parameterization can be assessed by comparing with the corre- 
sponding numbers of Ref. [E2|. For example, f w NN — 1-0 and f n NA = 1-8 in Ref. P2"| ; the 
effect of the parameterization is therefore of the order of 20%. 

Next, we calculated the full overlap function for p ^ 0. In Figure [| we plot the func- 
tion u(p) = O(p)/O(0) for the parameters given above. It is instructive to compare the 
momentum dependence of this form factor with the one given by the CBM |6,23[]: 



«(P) = 3^fl". ( 71 ) 

where jx is the spherical Bessel function and R is the radius of the underlying MIT bag. The 
solid line is our result and the dashed one is the CBM result for R = 1 fm. The faster falloff 
of our result is clearly a consequence of our Gaussian ansatz. As we will discuss soon, this 
rapid falloff will have the consequence of giving a smaller value of the self-energy correction 
to the nucleon mass, as compared to the corrections obtained with the CBM. 

VI. SELF-ENERGY CORRECTIONS TO THE NUCLEON AND A(1232) MASSES 



In this Section we present numerical results for the pionic self-energy corrections to the 
nucleon and A(1232) masses and discuss the LNA contribution to the masses. We start by 
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rewriting the vertex function in a manner to make clear the problem with the pion dispersion 
relation. The pion energy is given, for low k, as [ |14| . [25 



E 2 1 (k)=m 2 7r + k 2 



As) 

U (72) 



(t) 



where 



2 2m q (<ip^) 7 f dq . 

m - = / W \2 ' (TO = - 6 / 77^ sm( m) • ( 73 ) 



The point is that for an arbitrary quark-quark interaction one obtains in general two different 
values for the pion decay constant, /w and (the explicit calculations can be found in 



Refs. |14y25]), depending on how one defines the decay constant. When using the time 
component of the axial current, one obtains and when using the space component one 
obtains fjf\ However, as suggested in Ref. ||14|| , and explicitly demonstrated in Ref. |25 



this problem can be cured by adding a transverse gluon interaction. Therefore, to illustrate 
the point of obtaining the correct LNA term from Eq. fl52"|) with composite pions, we use 
the correct pion dispersion relation and assume fjp = = f w and denote Ei(k) = u(k). 
The normalization of the pion Salpeter amplitude, Eq. (0), can be rewritten as 

r rik 9 

A/» = MP) J 7^3 sin^(g) h(q) = -u>(p)f 2 . (74) 

That is, a 2 from Eq. ( |60"D is 2/3 f 2 . We next extract from the vertex function (we concentrate 
on the NN form factor) this normalization in the following way 

Wh N (p) = 1 G ^ P) Tj,<r N .p. (75) 
/2a; (p) 2 U 



The relation of the function G^ik 2 ) to the overlap function 0(p) can be trivially obtained 
by comparing with Eq. flST]) . 

Inserting Eq. ( [T5|) in the expression for the N and A self-energies, Eqs. (|52|) and (|53"D, 
and after performing rather straightforward spin-isospin algebra one obtains 



M N = i< 0) 



2 [°° p*u 2 { P ) 32 2 foo P 4 U 2 { P ) 

/o io dp ^W~^ U Jo u(p) [AM + u(p)) ' (76) 
MA = Mi + 25 /o io dp M\AM- M ] - f °Jo *"^P (77) 
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where 

AM = M A - M N , u(p) = yp 2 + ml, (78) 

and 

u(v) - f 2 - I?? IkK (79) 

U[P) ~ 0(0)' /o ~ m 2 4vr • l?9j 

Note that in principle we have different spatial dependencies for the NN, NA, • • • vertices, 
but for simplicity we have written them here as being equal. A schematic representation of 
Eqs. ([76]) and ( [77]) is presented in Fig. |3|. It is important to note that these equations are not 
the ones one would obtain by simple perturbation theory; they are actually nonperturbative, 
because of the dependence on AM = Ma — M^ on the r.h.s. 

It is easy now to obtain the LNA contributions to the masses ||. For the nucleon, the 
LNA contribution comes from the first term in Eq. ([7(]) by performing the integral. The 
integral can be done by transforming it into a contour integral and making use of Cauchy's 
theorem. The result is Eq. (0). For the A, the LNA contribution follows in a similar way 
from the last term in Eq. flT7|). 

To conclude, we discuss numerical results for the pionic corrections. Initially we solve 
variationally the bare nucleon case. As discussed above, using K$ = 247 MeV, we obtain 
for the variational size parameter the value = 1.2K Q . We also use here = a A . This 
leads to the following values for the bare A" and A masses 

il4 0) = 1 174 MeV , Mf = 1373 MeV . (80) 

The difference between the masses, of the order of 200 MeV, comes from the hyperfine 
splitting induced by the confining interaction. Given these values, we solve the two self- 
consistent equations given in Eqs. flTB|) and ([T7|). They are solved by iteration. We obtain 
for the masses 

M N = 1125 MeV , M A = 1342 MeV . (81) 
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Comparing with the values above, we see that the pionic effect is relatively small, as it 
should be, and of the order of 50 MeV for the N and 30 MeV for the A. The pionic effect is 
smaller for the A, as one expects from spin-isospin considerations B. The results obtained 



with the CBM for a R = 1 fm are a bit larger |23] . The difference can be traced to the rapid 
falloff of the form factor in our model. 

We certainly do not expect these numbers to be definitive. Once more realistic micro- 
scopic quark interactions and ansatze for the baryon wave function are used, they might 
be improved. However, independently of the microscopic model, our scheme is general and 
able to incorporate such interactions and new baryon amplitudes. I would be of particu- 
lar interest to have the numbers for a linear confining interaction with short range gluonic 
interactions that respect asymptotic freedom. 



VII. CONCLUSIONS AND FUTURE PERSPECTIVES 

We developed a calculational scheme to calculate chiral loop corrections to properties 
of composite baryons with composite pions. The composite baryons and pions are bound 
states derived from a microscopic chiral quark model inspired in Coulomb gauge QCD and 
provides a generalization of the Nambu-Jona-Lasinio model to include confinement and 
asymptotic freedom. An effective chiral hadronic model is constructed by projecting the 
microscopic quark Hamiltonian onto a Fock-space basis of single composite hadronic states. 
The composite pions and baryons are obtained from the same microscopic Hamiltonian that 
describes the chiral vacuum condensate. The projection of the quark Hamiltonian onto the 
pion states is nontrivial because of the two-component nature of the Salpeter amplitude 
of the pion. As explained before, the two components correspond to positive and negative 
energies which complicates the Fock-space representation of the pion state. The projection 
is made possible by rephrasing the formalism of the Salpeter equation in terms of the RPA 
equations. 

The development of models and calculational methods of the sort described in the present 



21 



paper are relevant in the context of a phenomenological understanding of nonperturbative 
phenomena of strong QCD like confinement and dynamical chiral symmetry breaking. Even- 
tually full lattice QCD simulations aimed at studying hadronic structure will be available 
and phenomenological models will play a central role in the interpretation of the data gener- 
ated. The developments of the present paper are of particular interest for the first-principle 
developments based on the QCD Hamiltonian, such as the nonperturbative renormalization 
program for the QCD Hamiltonian [|17| and Hamiltonian lattice QCD fl8 |. We intend to 
implement the technique developed here to such first-principle QCD calculations. 

We illustrated the applicability of the formalism with a numerical calculation using a 
simple microscopic interaction, namely a confining harmonic potential, and a simple Gaus- 
sian ansatz for the baryon amplitude. This very same S"-wave interaction has been used in a 
variety of earlier calculations, such as meson and baryon spectroscopy and S-w&ve nucleon- 
nucleon interaction. Numerical results were obtained here for the pion-nucleon form factor 
and for the pionic self-energy corrections to the nucleon and A(1232) masses in the non- 
perturbative one-loop approximation. Despite the simplicity of the interaction, the results 
obtained are very reasonable. 

For the future, the most pressing development would be to use a microscopic interaction 
that is consistent with asymptotic freedom and describes confinement by a linear potential. 
The calculation of the pion wave function beyond lowest order in momentum must be im- 
plemented and the variational ansatz for the baryon amplitude must be improved. A more 
ambitious development would be to include explicit gluonic degrees of freedom. In this case 
renormalization issues will show up and the new techniques such as discussed in Ref. Jl7| will 
certainly be useful. Another very interesting direction would be to employ the techniques 
developed here in Hamiltonian lattice QCD. 
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FIG. 1. Graphical representation of the functions Wf, i = 1, 2, 3. 
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FIG. 2. The function u(p). The solid line is the form factor obtained with the baryon amplitude 



of Eq. ( |57| ) and pion Salpeter amplitudes of Eqs. flog ) and (59). The dashed line is the CMB form 
factor of Eq. (|7l|) for R = 1 fm. 
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